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ABSTRACT. Self-assembly is a fundamental process by which supramolecular species form spon- 
taneously from their components. This process is ubiquitous throughout the life chemistry and is 
central to biological information processing. Algorithms for solving many mathematical and com- 
putational problems via tile self assembly has been proposed by many researchers in the last decade. 
In particular tile set for doing basic arithmetic of two inputs have been given. In this work we give 
tile set for doing basic arithmetic (addition, subtraction, multiplication) of n inputs and subsequently 
computing its modulo. We also present a tile set for primality testing. Finally we present a software 
'xtilemod' for doing modular arithmetic. This simplifies the task of creating the input files to xgrow 
simulator for doing basic (addition, subtraction, multiplication and division) as well as modular 
arithmetic of n inputs. Similar software for creating tile set for primality testing is also given. 

Keywords: DNA self-assembly, error-correction, algorithmic self-assembly, Wang tile, DNA comput- 
ing, Xgrow, XTile 

1 Introduction 

Self-assembly is a natural phenomenon observed at many places in nature such as formation 
of galaxies, formation of coral reefs, crystal growth etc. In 1996 Erik Winfree of California 
Institute of Technology, showed that it can be used to perform nano-scale computations 
[H5ll . This paved way for the birth of algorithmic self-assembly utilizing knowledge of three 
fields - DNA Nanotechnology [12J (due to the pioneering work of Ned Seeman in 1980s), 
DNA Computing (due to the pioneering work of L. Adleman in 1994) and Tiling Theory 
fl3~| (due to pioneering work of H. Wang who showed that zigsaw shaped colored tiles can 
simulate universal Turing machine). Winfree formulated the idea of molecular Wang tile 
using all this and showed that it can simulate universal Turing machine using abstract Tile 
Assembly Model (aTAM) ffl8l l8|. Algorithms for solving many mathematical and compu- 
tational problems via tile self assembly has been proposed by many researchers in the last 
decade (6J. |3 |4j, [7|. The addition of two numbers with Wang tiles is given in the book Q. 
In 2006 Brun fH 13 gave algorithm for addition and multiplication of two numbers with 
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aTAM and in 2008 Zhang and his coworkers [19J gave algorithm for subtraction and di- 
vision of two numbers in aTAM. In this work, we extend the results and provide tile set 
for addition and multiplication of n numbers together with a tile set for arithmetic expres- 
sion involving addition and subtraction. This is further extended to a tile set for modular 
arithmetic expression. Winfree's tile assembly model can be simulated using a program de- 
veloped by him called Xgrow fl7|. The input to Xgrow is a .tiles file. Many such standard 
.tiles files are distributed with the Xgrow package. In this paper, we also present a software 
package 'Xtilemod' for modular airthmetic that can be used to create such input files for 
the Xgrow simulator of Winfree by providing the basic arithmetic expressions. For further 
details on algorithmic self-assembly the reader is referred to excellent papers (161 ITT1 [10| 
and thesis |!T4ll5|. 

This paper is organized as follows. In Section 2, we give a brief overview of the aTAM. 
Section 3 provides the algorithms for addition and multiplication of n numbers via aTAM. 
Algorithm for arithmetic expression of n numbers involving addition and subtraction is 
described in Section 4 and algorithm for modular arithmetic expression of n numbers in- 
volving addition and subtraction is given in Section 5. Section 6 gives tile set for primality 
testing. In Section 7 we describe a software for modular arithmetic expressions which can 
produce the input to xgow simulator. Finally Section 8 concludes the paper with some gen- 
eral remarks. 

2 Background 

Self assembly of DNA molecules can be approximated by mathematical models. Several 
such models have been given by many authors. In [18|. Winfree gave a model called aTAM. 
The basic building block of a tile assembly model are square tiles. Each tile has 4 edges viz 
north, east, south and west with glues assigned to each edge. Each glue has some strength 
(usually an integer). Assembly starts with some seed tiles and frame tiles giving the initial 
input of the aTAM and a tile can attach to another tile if (1) the edges have matching glues 
and (2) the combined glue strength of both the tiles is greater than or equal to system tem- 
perature (again an integer, usually denoted by T and we take T = 2). Tiles can not be rotated. 
Under these settings system will grow from some initial tile set configuration to form a pat- 
tern or depending upon problem input it solves a problem. In Figure [lj a tile attaches to a 
seed configuration of 3 tiles. East glue a\ of configuration tile and west glue of outside tile 
is of strength 2, all other tiles have glues of strength 1 except north of configuration tile and 
south of each tile in seed configuration that has strength 0. System temperature is T = 2. It 
is shown by Winfree that such system is Turing universal. For more formal description the 
reader is referred to |3j[TJ[2]|. Recently basic arithmetic of two inputs with DNA Self Assem- 
bly have been studied by many authors OJ |2j [19). In particular, two algorithms have been 
proposed by Brun (TJ |2j for adding two numbers. Algorithm using L-configuration uses 
16 computational tiles for addition whereas there is another economical tile set for adding 
two numbers which takes just 8 computational tiles. In [2], Brun also gave algorithm for 
multiplication of two numbers with 28 computational tiles. Algorithm for subtraction and 
division are given in Ifl9l. In the next sections we will assume that the reader is familiar with 
the notations and definitions of |2lH9ll. 
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Figure 1: abstract Tile Assembly Model (aTAM) : A Tile attaches to a Seed Configuration of 
3 Tiles 

3 Basic Arithmetic of n inputs with DNA Self Assembly 

3.1 Addition of n Numbers with 8 Computational Tiles 

Suppose we want to add n positive integers {fli,fl2, ■ . - ,an\ai £ Z + }. We can write their 
binary representation as cij = Y^gLo a 0^ \ where nij = \aj\ is the size of the integer fly and 
ccg £ Z2. If m = max{nij\l < j < n} then we can pad extra zeros to write fly = X^5Lo ^e^' . 
The next two Theorems gives the tile set for computing Yd=i a i- 

Theorem 1. Let 

Ex = {s,£,0,l,00,01,10,ll}U{#z|l < i < (m + n)(n-2)} 
U{$z'|l < i < max(m + n,2n -2)} 

be the set of glues and let T\ be a set of tiles (different tiles types of T\ are given in Table^) 
over Ei as described in Figure^then T\ computes the function £f =1 fl,-. 

The logic of the system is identical to a series of one-bit full adders. Each solution tile 
takes a bit from each of the inputs on the south side and a carry bit of the previous solution 
from the east side, and outputs the next carry bit on the west side and the sum on the 
north side. After computing the sum of two numbers, we need to combine this sum with 
next input bitwise. This process continues till all the integer inputs are added. To handle 
overflow we make size of each input fl ; - as n + m — 1 bits. We will illustrate it by an example. 
Figure |3] shows a sample addition of 12 + 6 + 2 + 4 = 24 (m = 4, n = 4). Input 12 = 1100 2 
(size 4) and 6 = IIO2 (size 3) are given at the bottommost (first) row by adding extra zeros 
(to adjust overflow we make size of each n + m — 1 = 7 bits) and reading 12 from right side 
and 6 from left side. The sum is computed (12 + 6 = 18 = IOOIO2) in second row and then 
passed to the 3 rd row from left side, next input 2 = IO2 is given from right side of 3 rd row 
and the sum is computed in 4 row (12 + 6 + 2 = 20 = IOIOO2) which is then passed to the 
5 th row from left side, next input 4 = IOO2 is given from right side of 5 th row and the final 
sum 12 + 6 + 2 + 4 = 24 = OOIIOOO2 is computed in 6 th row. Last row is added to stop the 
computation. Note that the example uses 5 left frame tiles, 4 corner tiles, 5 right frame tiles, 



4 

7 top frame tiles, 21 input tiles and 14 computational tiles (overall) matching the data given 
in Tabled] 

Remarks: We have used 8 computational tiles in the algorithms, however we can optimize 
it a bit since the lower portion (10 or 01) of computational tiles gives same effect in com- 
puting the sum. Hence while implementing it in "xtilemod" (see Section [7]| we have used 6 
computational tiles. 

Table 1: Total Tiles used for Computing Yli=i a i using 8 Computational Tiles 



Basic Tiles 


Tiles Types 


Overall Tiles 


Left-Frame 


n-1 


2n-3 


Corner 


4 


4 


Right-Frame 


2n-3 


2n-3 


Top Frame 


2 


n+m-1 


Input 


(2n-3)(n+m-l) 


(n+m-l)(n-l) 


Computational 


8 


(n+m-1 )(n-2) 



3.2 Addition of n Numbers with L type Tiles 

Theorem 2. Let E 2 = {0,l,s,ss} U {$z'|l < i < n + tn- 1} U {#c|0 < c < n} and T 2 be a 
set of tiles (different tiles types of T2 are given in Table^) over E2 as described in Figure^ 
Then T2 computes the function Y%=\ a i- 

The logic of this system is similar to Theorem 2.2 of O. The tile set of addition of n 
numbers using L type tiles of [2 J have been implemented in xtilemod (see Section [7]). 

Table 2: Total Tiles used for Computing Ya=i a i using L-type Computational Tiles 



Basic Tiles 


Tiles Types 


Remark 


Left-Frame 






Right-Frame 


EL2 m k 




Top Frame 


2+ (n + m-1) 




Computational 




Bruns' L-Type Tiles 



3.3 Multiplication of n Numbers with Tiles 

Theorem 3. LefZ 3 = {0, 1,00, 01, 10, 11, 20, 21, #} U {#/l < i < nm-2}U{$; : nm < j < 
nm + (n — 1) + Y%=2 m k} ^3 ^ e a se ^ °^ ^ es (different tiles types of T3 are given in Table 
[3] ) over E3 as described in Figure^ Then T3 computes the function nf=i a i- 

Note that the size of nf=i a i 1S n + m bits, where m = max{m ; |l < j < n}. WLOG, 
Suppose fli is of max size m. We put a\ as an input at the bottommost horizontal base tiles 
and all other inputs are given in the vertical frame (right column) followed by a separator 
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Figure 2: Til e Set for Com puting £" =1 flj using 8 Computational Tiles. Note that in Input 
Tiles (Figure 2(g) and 2(h) i each oc ] e G Z2 also a, b, c G {0, 1}. 




n+m-1 



Figure 3: Tiles Growth for Expression 12 + 6 + 2 + 4 = 24 using 8 Tile Types Configuration 



tile. The idea of multiplication of n input integers is a simple extension of the Theorem 2.4 
of dJ. The computational tiles also remain the same. The only difference is that n — 2 inputs 
are given in the vertical frame separated by one separator tile and subsequently one row. 
Figure [6] shows the product of 5 x 4 x 3 = 60. Note that the input 5 = 1OI2 is given at the 
bottom row, input 4 = 1OO2 is given as the first vertical input and 3 = II2 is given as second 
vertical input. The output is given at the top row as 60 = IIIIOO2. 



Table 3: Total Tiles used for Computing Yl?=i a i 



Basic Tiles 


Tiles Types 


Remark 


Horizontal-Frame 


2(nm - 1) 




Seed 


1 




Vertical-Frame 


2n + rrij 




Top Frame 


6 




Computational (including a copy tile) 


9 





4 Arithmetic Expression of n Numbers with Addition and Sub- 
traction: E?=lft«i/ & e { + 

Theorem 4. LerZ 4 = {0, 1, 00, 01, 10, 11, c 0/ c_i} U {d,e,f,g r h,l} be the set of glues and let 
T4 be a set of tiles (different tiles types of are given in TaWe^) over E4 as described in 
Figure^ Then T4 computes the function Yd=\ fii a u Pi £ {+1/ — !}• 
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(f) Corner Tiles 



(g) Example showing addition of 6 + 4 + 3 + 5 = 18 



Figure 4: Tile Set for Computing Y%=\ a i using L Type Computational Tiles. Note a,b,c G Z2. 
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(c) Vertical Frame Tiles (d) Seed and Result Tiles 

Figure 5: Tiles set for Computing Ylf=i a i- Note a, b, c £ Zz- 




Figure 6: Growth of Tiles for Computing the Product of 5 x 4 x 3 = 60 



Preprint 2012 



9 



The logic of the system is identical to a series of one-bit full adders or one bit subtractor 
depending on the operation which is to be performed. Each solution tile takes in a bit from 
each of the inputs on the south side and a carry bit from the previous solution tile on the 
east side, and outputs the next carry bit on the west side and the sum/ subtracted value 
on the north side. In the subtraction the carry bit which is generated on the west indicates 
the requirement of 1 from the next left tile attached to it. For separating the operation of 
addition and subtraction + or - sign is attached in all the computational tiles on the east and 
the west side along with the computational bit. Due to this addition and subtraction can 
be done by just supplying the sign to computation from the right frame tile. We need to 
combine this result with next input bitwise. This process continues till all the integer inputs 
are operated. After operating it if output is positive then the same is displayed on the top 
row(which is the final output) but if output is negative then the output is converted into 
unsigned binary bit and this will be displayed on the top row(which is the final output) 
from the 2s complement. This is done because the normal output which was generated 
for negative value is in the 2s complement form. Also a negative sign is displayed with 
the negative result to tell the user that final output is negative. We will illustrate it by an 
example. Figure [8] shows a sample addition of 6 — 12 + 4 — 2 = —4. Input 6 = IIO2 and 
12 = IIOO2 are given at the bottommost (first) row by adding extra zeros and reading 6 
from right side and 12 from left side. The subtraction is computed (6 — 12 = —6 = IIOIO2) 
in second row and then passed to the 3 rd row from left side, next input 4 = IOO2 is given 
from right side of 3 rd row and the sum is computed in 4 row (6 — 12 + 4 = —2 = IIIIO2) 
which is then passed to the 5 th row from left side, next input 2 = IO2 is given from right side 
of 5 th row and the final result 6 — 12 + 4 — 2 = —4 = IIIIIOO2 is computed in 6 th row. Now 
as output is negative we have to convert it back into normal unsigned binary expression 
from 2s complement with a negative sign in front of it. To this the row 7 th ,8 th ,9 th are used. 
Also the final result — OOOOIOO2 is printed on the 9 row. 



Table 4: Total Tiles used for Computing YH=\ fii a i> Pi e {+1/ — 1} 
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2n+5 




Corner 


6 
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2n-l 




Top Frame 


6 




Input 


(2n-3)(n+m-l) 




Computational (+) 


8 




Computational (— ) 


8 




Other Computational 


4 





10 



a©b0c- 



= aeb©c + 



a b 



J 



1 11 ir'i 



I 1 I 1 I 1 1 I 1 1 

Li oJ Li cJ Li J L J 



(a) Typical Computational Tiles (b) All 8 Computational Tiles for Addition 
for Addition, where 

cout = (a A b) V (b A c) V (c A a) 

sub 
sub 



n 



u l'ii 1 r ^ 

l — 1 I II 1 II 1 I 

I , oJ L, oJ L, J L ,J 

(c) Typical Computational Tiles for Subtrac- (d) All 8 Computational Tiles for Subtraction 
tion (Refer Table [6} 




j a i-i a J" 

v>. i 2 

+ a i-i a j-i ¥V 



(e) Input Tiles for Bottom Row for Two Inputs 

l<i<m + n — 1 



S2i =, 



p 



ER 

$2i- 



,= S2'i+1= 


r — 


r = 


£ SR 


CL 




1= S2i = 


l=S2{n-l) = 





(g) Right Frame Tiles (ER : 1 < i < n — 1) and 

(SR : 1 < i < n -2) and A = 1 + (/ - l)(n + m), 0' 6 {+,-} 

depending on position i 



| "'■-I — | | 1 1 

LoJLj 

(f) Intermediate Input Tiles (1 < z < n + m - 
1),3 < < n and A = i + (j — 3)(n + m) 



Q_ 

1-1 

n 

(h) Left Frame Tiles: 1< i < n - 2), s e {CO, C_l} 



m r$1 " 



o 


r ' -i 


u 


EL 


SL ! 




L s j 


_ , J 






- co n 


u 


EL ° 


^^■ + 


EL " 


L. J 


— C_l — 1 


— CO — 



L L J 



JL.JL JL.J 

(i) Top Frame Tiles 

mm 

a. I a H c. I 

L ,JL, 



JL J 



e = d 

(j) Corner Tiles 




a 



7- t;i^ c^t- t^v r r 



(k) Other Computational Tiles 



Preprint 2012 




— C_l ■ 
EL 

|C_1 | 

SL 
|C_1 | 

EL 

| C_l | 

SL 

' C_l | 
E_ 
CO 

S. 







a 
i 
l 

- i 
l 

- 01 
01 

- 1 
1 

- 01 
01 

I 1 
1 

■ 00 



+ ■ + 



1 

I 1 
1 

- Dl 
01 

- 1 
1 

- Dl 






1 
1 

1 

- Dl 
01 

| 1 
1 

- 01 





I 


I 1 
1 

1 

I 01 
01 

I 1 
1 

- 01 



3 oi I oi 5 oi 



1 

10 ■ 



1 

I 



1 I 
1 

1 I 
1 

Dl - 
01 

I I 
1 

10 ■ 
10 

o I 



II ■ 



°+° + 

Hal 

+ 

3 

+ 

00 — 1= 

8 OO 5 

1 ol 

+ 



SB 

:• er 

MJ — 1= 55 

00 3 SR 
,;, _f= S4- 

.:■ 1 er 

S3 



01 5 00 

1 





00 ■ 



SR 
52 : 
ER 

51 : 



oo S oo Sr oo 



11 a oi B oo 



Figure 8: Example showing 6 — 12 + 4 — 2 = —4 

5 Modular Arithmetic Expression of n Numbers with Addition 
and Subtraction: EJLi fa a u (mod t) e {+1, -1} 

The tile set for computing Y%=\ ¥>i a i' (mod t) jS, S {+1/ —1} is just an extended version 
of the Yd=i C>i a if Pi S {+1, —1} tile set (See Theorem[4]). We add tile set for division (from 
Lemma 2 of HH) to the output of the £f =1 jSjflj, /3, G {+1, —1} tile set after ignoring sign bit. 
The final remainder of the output comes on the top most row. An example for 6 — 12 + 4 — 2 
(mod 3) = 1 is given in Figure|9] An implementation of the tile set is given in "xtilemod". 



6 Primality Testing with Tiles 

Theorem 5. Let Z 5 = {0, 1, 000, 001, 010, 011, 100, 101, 110, 111} U {= 0, = 1, = 2, 0-, 1-, = 
, >, <, SS, RR,JJ, TT, WW, ZZ, KK, YY, XX, Y, K, X} be the set of glues and let T 6 be a set of 
tiles (different tiles types ofT$ are given in Table^)over E5 as described in Figure\i0\ Then 
T5 determines of primality of any input integer n. 

The logic of primality testing tile set is based on the Algorithm [l] We will illustrate it 
by an example. Figure 11 shows a sample primality testing for n = 5 = IOI2. This binary 
input is given at the bottommost (first) row. These bits are copied to left side of the centre of 
each tile of 3 bits of second row. Second row also computes k = |§ J and stores it in the right 
side of the centre of 3 bits of each tile. The middle bit out of these 3 bits on each tile centre 
in the second row is intermediate computation bit I (which is initially the bit of n). In fact 
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Figure 9: Example showing 6 — 12 + 4 — 2 = —4 modulo 3 
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our input n always remain on left side in each row at the centre of each tile and the value of 
k remains same till I — k > is true otherwise it changes to k — 1 and I becomes n. Third 
row checks if A: > 1, which is true in this case so we copy the content (all 3 bits) of centre 
of each tile to 4 th row (if k is not > 1 it quits and displays that it is not a prime). Next row 
5 th checks if the content I corresponding to middle bit is > k (which is true) so I = I — k is 
done in 6 th row. Next row further checks if I — k > 0, which is true so we repeat previous 
step and I = I — k is done in 8 th row. Now 9 th row again checks if I — k > which is not 
true this time so in 10 th row we do k = k — 1 and I = n. In 11 th row again we check if k > 1 
(just like 3 rd row) and we find that k = 1 so 12 th row gives the final output as 5 is prime. 

Algorithm 1 Primality Testing Algorithm 
Require: n 
Ensure: k = [j\ 
while k ^ 1 do 

if n (mod k) = then 
Not Prime 
exit 
else 

k:=k-l 
end if 
end while 
Prime 



Table 5: Total Tiles used for Primality Testing 
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7 Xtilemod: A Software Package for Modular Arithmetic Expres- 
sions 

Xtilemod is an modular arithmetic software package implemented in java which generates 
.tiles file according to options selected at runtime. For example if we have to add two num- 
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(d) Right Frame Tiles 
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able 6: Subtraction Table Bitwise 
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;ure 11: Primality Testing Example for p = 5 
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bers using 8 tile type addition then we will select the option 1 which chooses the operation 
of two numbers and thereafter option 1 of 1 which specifies 8-tile type addition. Last step to 
get the resultant tile design is to enter the input of integers. The code will display the name 
of file where the resultant output is loaded. For example if we are adding 6 and 12 using 8 
tile type addition the output will be loaded in 

add_8_tile_12, 6 . tiles . 
The resultant .tile file can be simulated using xgrow. The display will be binary. Red color 
of the resultant tile indicates 1 and the white indicates 0. Resultant tiles are generally at 
north with a few exception. In case of the division tile set the xgrow compilation displays 
the quotient at the west and the remainder at the north. In 8 tile addition of two numbers 
the output is shown at the middle and the inputs are both in the north and south direc- 
tion. For more details the reader is referred to the user manual of the software available 
at |http : / / www.g uptalab .org / xtilemod/manual.pdf . The tool 'xtilemod' is available for 
download at http:/ / www.guptalab.org/xtilemod. 

8 Conclusion 

In this work we presented tile set for modular arithmetic expressions involving addition, 
subtraction, multiplication of n integers followed by a modulo operation and a software 
"xtilemod" for generating the corresponding .tiles files. A tile set for primality testing is 
also presented together with a java program "prime.java" for generating .tiles file. Due to 
space constraints we have omitted the proofs but the key ideas of each result are illustrated 
by examples, the detailed proof of the results will be reported in an extended version of the 
paper elsewhere. 
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